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We report a new measurement of the time-dependent CP-violating parameters in B° —* 7r + 7r~ 
decays with 535 x 10 6 BB pairs collected with the Belle detector at the KEKB asymmetric-energy 
e + e~ collider operating at the T(45) resonance. We find 1464 ± 65 B° — > it + tt~ events and 
measure the CP-violating parameters S^-k = —0.61 ± O.lO(stat) ± 0.04(syst) and Amr = +0.55 ± 
0.08(stat) ± 0.05(syst). We observe large direct CP-violation with a significance greater than 5 
standard deviations for any Snn value. Using isospin relations, we measure the Cabibbo-Kobayashi- 
Maskawa quark-mixing matrix angle </>2 = (97 ± 11)° for the solution consistent with the standard 
model and exclude the range 11° < 4>2 < 79° at the 95% confidence level. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 



In the standard model (SM) framework, CP viola- 
tion is attributed to an irreducible complex phase in the 
Cabibbo-Kobayashi-Maskawa (CKM) weak-interaction 
quark-mixing matrix Q. In the decay chain of T(4S*) -> 
B°B°, one B° decays into tt + it~ at time t v „, while the 
other decays at time f ta g into a flavor specific state / tag - 
The time-dependent CP violation Q is given as 



>-|A*l/T-flO 



[1 + q ■ {<S W7r sin(AmrfA<) 



4r B o 

-A^ cos(Am d At)}] 



(1) 



where At = — ttag, t b o is the B lifetime, Am^ is the 
mass difference between the two B mass eigenstates Q 
and q = +1 (—1) when / tag = B°(B°). <S W7r and A*-* are 
the mixing-induced and direct CP-violating parameters, 
respectively. 

The CP- violating parameters have been measured by 
the Belle Q and BaBar [j| collaborations. Both exper- 
iments obtained consistent results for S^ n . In contrast, 
BaBar measured an A^ value consistent with zero, while 
Belle found evidence for large direct CP violation with a 
significance of four standard deviations (a) using a data 
sample containing 275 x 10 6 BB pairs. Here we report 
a new measurement with a large data sample (535 x 10 6 
BB pairs), and improvements to the analysis method 
that increase its sensitivity. We confirm our earlier re- 

| 



decays at the 5.5<r level; the disagreement with the BaBar 
A W ir measurement [H remains. This result rules out or 
strongly constrains superweak models 0, extensions of 
the SM in which all CP violation occurs through AB = 2 
processes. 

One of the CKM angles, 4>2 0, can be measured using 
(Stttt = yl — A^n sin(2ri)2 + «0, where k is determined 
using isospin relations [9( . This angle has been measured 
using not only B — > irir decays but also B — > pit and pp 
decays |Io| ; all the measurements give consistent results, 
and the combined fa value, together with measurements 
of other CKM angles and sides, is consistent with the 

We combine our 
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unitarity of the matrix [11 . 
and A-xir measurements with the world average (W.A.) 
values of other quantities to obtain a new constraint from 
B — » 7T7T on fa. Multiple solutions are found; for the 
solution consistent with other CKM measurements in the 
context of the SM, the constraint is more restrictive than 
those obtained from other B decay modes. 

The data sample used in this analysis was col- 
lected with the Belle detector [Hj at the KEKB e+e" 
asymmetric-energy (3.5 on 8 GeV) collider Q operating 
at the T(4S') resonance produced with a Lorentz boost 
factor of (3j = 0.425 nearly along the electron beam di- 
rection (z axis). Since the two B mesons are produced 
approximately at rest in the T(45) center-of-mass sys- 
tem (CMS), the decay time difference At is determined 
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from the distance between the two B meson decay ver- 
tices along the z-direction (Az): At = Az/c(3j. In the 
Belle detector, a silicon vertex detector and a 50-layer 
central drift chamber (CDC) are used for charged parti- 
cle tracking, and an array of aerogel threshold Chercnkov 
counters as well as the dE/dx measurements in the CDC 
provide the particle identification (PID) information to 
distinguish charged pions and kaons. The devices are 
placed inside a superconducting solenoid coil providing a 
1.5 T magnetic field. 

We employ the event selection of Rcf. Q except for the 
PID requirement, which is removed. This increases the 
signal detection efficiency by 23 %. The PID information 
is instead used in a likelihood fit in this analysis, improv- 
ing the measurement errors for the CP- violating parame- 
ters by about 10 % compared with the previous analysis. 
We reconstruct B° — > 7r+7r~ candidates using oppositely 
charged track pairs. We select B meson candidates using 
the energy difference AE = E B —E^ cam and the beam en- 
ergy constrained mass M hc = yJ(E£ cam ) 2 - {p* B ) 2 , where 
E'beam is the CMS beam-energy, and E* B and p* B are 
the CMS energy and momentum of the B candidate. 
We define the signal box as 5.271 GeV/c 2 < M bc < 
5.287 GcV/c 2 and \AE\ < 0.064 GcV. 

The standard Belle algorithm in Ref. [IH identifies the 
flavor of /tag using properties of its decay products, and 
provides two variables: q defined in Eq. |T]) and 7'. The 
parameter r ranges from r = indicating no flavor dis- 
crimination to r = 1 indicating unambiguous flavor as- 
signment. The candidate events are categorized into six 
r intervals (I = 1,6). The wrong tag fraction in each I 
bin, wi, and the differences between B° and B° decays, 
Awi , are determined using data [l5| . 

To discriminate the continuum background (e + e~ — > 
qq, q = u,d,s,c), we form signal (qq background) like- 
lihood functions, Cs(bg), from features of the event 
topology, and require r-dependent thresholds of 1Z = 
£-s/(£s + Cbg) for the candidates, as the separation 
between signal and qq background depends on r. The 
thresholds are determined to be 0.50, 0.45, 0.45, 0.45, 
0.45, and 0.20 for each I bin by optimizing the expected 
sensitivity using signal Monte Carlo (MC) events and 
events in the sideband region 5.20 GeV/c 2 < Mb c < 
5.26 GeV/c 2 or +0.1 GeV < AE < +0.5 GeV. We fur- 
ther divide the data sample into two categories having 1Z 
above or below 0.85 to take into account the correlation 
between the AE shape of qq background and 1Z. We thus 
have 12 distinct bins of 1Z vs. r; these bins are labeled 
I = 1, 6 (I = 7, 12) for the six r intervals with 1Z > 0.85 
(1Z < 0.85). 

We extract signal candidates in the global area Mbc > 
5.20 GeV/c 2 and -0.3 GeV< AE < +0.5 GeV by apply- 
ing the above requirements and the vertex reconstruc- 
tion algorithm used in Ref. The selected candi- 
dates are not only B° —> tt + it~ signal events but also 
include B° — > K + ir~, qq and three-body B decay back- 



grounds. We estimate the signal and background yields 
with an unbinncd extended maximum likelihood fit, mak- 
ing use of AE, Mb c and the kaon identification probabil- 
ity x± = Cr± /{Ck± + £?r±) for the positively and neg- 
atively charged tracks of B° — ► tt + tt~ candidates, where 
C 7I ± and Ck± are the likelihood values for the pion and 
kaon hypotheses. 

We use a sum of two bifurcated Gaussians and a sin- 
gle Gaussian to model the AE and Md C shapes, respec- 
tively, for both B° — > ir + ir~ and K + tt~. The proba- 
bility density functions (PDF) as a function of x± for 
the signal and B Q — > K + tt~ decays are obtained from 
a large data sample of inclusive D* + — > D a ir + ,D — > 
K~ir + decays. The yields of B° — > K + ir~~ events are 
parameterized as n K ±^ = riKw(l T AftT7r)/2, where 
A K * = -0.113 ± 0.020 3 is the direct CP asymme- 
try in B° — > K + tt~ decays. We fix the Ak-k value and 
float the B° — > K + ir~ yield in the fit. For the signal 
and B° — » K + tt~ events, we use MC-detcrmincd event 
fractions in each JZ-r bin. 

The qq background shapes in AE and M^ c are de- 
scribed by a second-order polynomial and an ARGUS 
function [l7|, respectively. We use gr-dependent two- 
dimensional (x+ ,X-) PDFs for the qq background to take 
into account the correlation between positively and neg- 
atively charged tracks. These PDFs are determined from 
the sideband events. 

For the three-body B decays, we employ a smoothed 
two-dimensional histogram obtained from a large MC 
sample for the AE-My, c shape. We use the same x± 
PDFs as those of the signal and B° — > K + ir~ decays, 
with a Ai?-dependent kaon fraction determined from the 
MC sample. 

By fitting to the data in the global area we deter- 
mine the yields of the signal and background components. 
Interpolating to the signal box, we obtain 1464 ± 65 
7T+7T-, 4603 ± 105 K+tt- and 10764 ± 33 qq events, 
where the errors are statistical only. The contribution 
from three-body B decays is negligible in the signal 
box. From the signal yield and the detection efficiency 
(53.1 %), we estimate the measured branching fraction 
to be (5. 2 ± 0.2) x 10~ 6 , in agreement with the W.A. 
value [18(. Figure [T] shows the projection plots of AE, 
Mbc and x± for candidate events. 

To determine and A^, we apply an unbinned 
maximum likelihood fit to the A< distribution of the 
16831 candidates in the signal box. The signal distri- 
bution in Eq. ([1]) is modified to incorporate the effect of 
incorrect flavor assignment, using wi and Awi. This dis- 
tribution is then convolved with the proper time interval 
resolution function R s i g (At) The final signal PDF 

is given by P^(Ai) = (1 - U)P&(At) ® R sig (At) + 
foiPoi(At), where the outlier PDF V i(At) accommo- 
dates a small fraction f Q \ of events having large At values. 
The At distribution for B° -> K+ir' is T"^ ±7iT (At) = 
(l/4T B o)e~l At l/ T s°[l - qAwi + q(l - 2w t ) cos(Am d Ai)]; 
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1: (a) AE, (b) Mbc and (c) x± projection plots of the 
* 7r + 7r~ candidates having 1Z > 0.85 in the signal box of 
(a) M bc with x± < 0.4, (b) AE with x± < 0.4 and (c) M bc 
with GeV< AE < 0.02 GeV. Figure (c) is the sum of x+ 
and X- distributions. 



the corresponding PDF P K ± T (At) is constructed in the 
same manner as the signal PDF. The qq background 
distribution contains prompt and finite-lifetime compo- 
nents; it is convolved with a background resolution func- 
tion modeled as a sum of two Gaussians and combined 
with the outlier PDF to give the qq background PDF 
Pqq(At). All the parameters of P q q(At) are determined 
using sideband events. 

We define a likelihood value for the i-th event, which 
lies in the ^-th bin of 1Z vs. r: 



p=Y.<n w &)pr\^) 



(2) 



is the fraction of component k G 



Here 

{TT + Tr-,K + TT-J\-n + ,qq} in TZ-r bin l\ V% W \s) 



the event-by-event probability for component k as a 
function of s = (AE, M bc , x+, x_); and P^ (At) is the 
event-by-event probability for component k and flavor 
tag q as a function of At. In the fit, 5 W7r and are the 
only free parameters and are determined by maximizing 
the likelihood function C = \\ i Pi . 

The unbinned maximum likelihood fit yields = 
-0.61 ± O.lO(stat) ± 0.04(syst) and A^ = +0.55 ± 
0.08(stat) ±0.05(syst). The correlation between S^n and 
Atttt is p = +0.15. Figures [H (a) and (b) show the back- 
ground subtracted At distributions of the signal events 
with r > 0.5 for q = ±1 and the asymmetry Acp in each 
At bin, respectively, where A C p = (N + -N-)/(N + +NJ) 
and N + r_\ is the number of signal events with q = +1 
(— 1) obtained by a fit in each At bin. 

The main contributions to the systematic error are due 
to uncertainties in the vertex reconstruction (±0.03 for 
Stttt and ±0.01 for A nv ) and event fractions (±0.01 for 
iStttt and ±0.04 for A^); the latter includes a conser- 
vative uncertainty for the possible qq background flavor 
asymmetry of ±0.02. We include the effect of tag side in- 
terference [H on (±0.01) and A^ (±0.02). Other 
sources of systematic error are the uncertainties in the 
wrong tag fraction (±0.01 for both and A^), physics 
parameters (tbo, Amj and Ak-k) (< 0-01 for both 
and Ann), resolution function (±0.02 for both 5^ and 
Atttt), background At shape (< 0.01 for both <S ffff and 




FIG. 2: (a) At distributions of B° — > tt + -k~ signal events 
with r > 0.5 after background subtraction for q — +1 (solid) 
and q = — 1 (dashed), and (b) asymmetry Acp plot. The 
curves are projections of the fit result. The difference in the 
heights of the q = +1 and q = — 1 components in (a) is due 
to direct CP violation. 



Atttt), and fit bias (±0.01 for both 5 WT and A^). We 
add each contribution in quadrature to obtain the total 
systematic error. 

To validate our CP-violating parameter measurement, 
we check the measurement of A^ using a time- integrated 
fit, and obtain A^ = +0.56 ± 0.10, consistent with the 
time-dependent fit results. An unbinned extended max- 
imum likelihood fit to the q = +1 (q = —1) subset with 
TZ > 0.85 and r > 0.5 yields 280 ± 20 (169 ± 16) tt+tt" 
signal events, in agreement with the measured A^n value 
taking into account the dilution due to the wrong tag 
fractions and B a B a mixing. We also check the direct 
CP asymmetry in B° — > K + n~~ events by floating Ak-k 
in the time-dependent fit, and obtain a value consistent 
with the W.A. Q and the same p value with the nominal 
fit. The fit is applied to various data subsets: a subset 
containing events with positive (negative) AE in which 
the B — > K + 7r~ contamination is suppressed (enriched), 
where A™ = +0.60 ± 0.11 (+0.51 ± 0.12), events with 
TZ > 0.85 (7Z < 0.85) where the qq background fraction 
is suppressed (enriched), events with x± < 0.4 where the 
signal fraction is enhanced, and events in one of the six r 
bins having different wrong tag fractions. All fits to the 
subsets yield CP asymmetries consistent with the overall 
fit result. We also carry out a fit to the sideband events, 
and find no sizable asymmetry. 

We determine the statistical significance of our mea- 
surement using a frequentist approach [2l|, taking into 
account both statistical and systematic uncertainties. 
Figure [3] shows the resulting two-dimensional confidence 
regions in the S^n and Ann plane. The case of no direct 
CP violation, A^ir = 0, is ruled out at a confidence level 
(C.L.) of 1 — 4 x 10 -8 , equivalent to a 5.5c significance 
for one-dimensional Gaussian errors. We also observe 
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FIG. 3: Confidence regions for S^n and Aw The curves 
show the contour for 1 - C.L. = 3.17 x KT 1 (Iff), 4.55 x 10~ 2 
(2a), 2.70 x 10" 3 (3ff), 6.33 x 1(T 5 (4a), 5.73 x KT 7 (5ff) and 
1.97 x 10 -9 (6ff) from inside to outside. The point with error 
bars is the S^tt and An* measurement. 



mixing- induced CP violation with a significance greater 
than 5.3cr for any Ann value. 

To constrain (f> 2 , we use the isospin relations Q with 
our measured values of S^tt and A^n the W.A. branching 
ratios of B° -> tt+tt" (5.2 ± 0.2 in units of 10~ 6 ), 7r°7r° 
(1.31 ± 0.21) and B+ -> tt+tt (5.7 ± 0.4), and the W.A. 
direct CP-asymmetry for B a -> 7r°7r° (0.36 ± 0.33) [II]. 
The isospin relations imply that the branching ratios and 
CP-violating parameters can be expressed with six pa- 
rameters, one of which is (f> 2 . We construct a x 2 using 
the measured values and the six parameters. We fol- 
low the statistical method of Ref. [ll|. We minimize 
the X 2 (4>2) for each (j> 2 value by varying the remain- 
ing five parameters with constraints on the correspond- 
ing branching ratios and asymmetries and compute the 
C.L. from the cumulative distribution of the difference 
&x 2 = X 2 (02) — Xmin fo f one degree of freedom, where 
Xmin i s the minimum x 2 { ( t ) 2) value. Figure [4] shows the 
difference 1— C.L. plotted for a range of <f> 2 values. We 
find four different solutions consistent with our measure- 
ment. For the solution consistent with the expectation 
from other CKM measurements, (lOO^) [H|], we nnd 
<t> 2 = (97±11)°. We exclude the (j> 2 range 11° < (j) 2 < 79° 
at the 95% confidence level. 

In summary, using a data sample containing 
535xl0 6 PP pairs we measure the CP- violating param- 
eters in B a — > 7r + 7r~ decays: S^-n = — 0.61 ± O.lO(stat) ± 
0.04(syst) and A™ = +0.55±0.08(stat)±0.05(syst). We 
report the first observation of direct CP violation with 
5.5cr significance. Our results as well as the evidence for 
direct CP violation in B° — > K + ir~ decays [22[ rule out 
superweak models @- The measured 5 W7r and A^-n val- 
ues in this Letter are consistent with those reported in 
Ref. [ij], and supersede Belle's earlier evidence for direct 




30 60 90 120 150 180 
4> 2 (degrees) 

FIG. 4: Difference 1-C.L. plotted for a range of <j>2 values 
obtained with an isospin analysis using Belle measurements 
of 5-n-n- and A-n-w, and the W.A. values for the direct CP 
asymmetry in B° — » 7r°7r° decays and branching fractions 
for the B — * nn modes. The solid and dashed lines indicate 
C.L. =68.3% and 95%, respectively. 



CP violation. Among the four (j> 2 solutions, the ±ler 
range for the 4> 2 solution consistent with the SM is more 
restrictive than that from measurements of B — > pir and 
B — > pp decays fiol ] . 
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